A solvable model which has X(5) as a limiting symmetry and 
removes some inherent drawbacks 
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^ Abstract 

Solvable Hamiltonians for the (3 and 7 intrinsic shape coordinates are proposed. The eigenfunc- 
tions of the 7 Hamiltonian are spheroidal periodic functions, while the Hamiltonian for the /3 degree 
of freedom involves the Davidson's potential and admits eigenfunctions which can be expressed in 
terms of the generalized Legendre polynomials. The proposed model goes to ^(5) in the limit 
of |7|-small. Some drawbacks of the X(5) model, as are the eigenfunction periodicity and the 7 
Hamiltonian hermiticity, are absent in the present approach. Results of numerical applications to 
^^^Nd, ^^^Gd and ^^^Os are in good agreement to the experimental data. Comparison with -^^(5) 
calculations suggests that the present approach provides a quantitative better description of the 
data. This is especially true for the excitation energies in the gamma band. 



PACS numbers: : 21.10.Re, 23.20.Lv, 21.60. Ev 



I. INTRODUCTION 

Since the liquid drop model was developed the quadrupole shape coordinates were 
widely used by both phenomenological and microscopic formalisms to describe the basic 
properties of nuclear systems. Based on these coordinates, one defines quadrupole boson 
operators in terms of which model Hamiltonians and transition operators are defined. Since 
the original spherical harmonic liquid drop model was able to describe only a small amount of 
data for spherical nuclei, several improvements have been added. Thus, the Bohr-Mottelson 
model was generalized by Faessler and Greiner 2| in order to describe the small oscillations 
around a deformed shape which results in obtaining a flexible model, called vibration ro- 
tation model, suitable for the description of deformed nuclei. Later on this picture was 
extended by including anharmonicities as low order invariant polynomials in the quadrupole 
coordinates. With a suitable choice of the parameters involved in the model Hamiltonian the 
equipotential energy surface may exhibit several types of minima Q] like spherical, deformed 
prolate, deformed oblate, deformed triaxial, etc. To each equilibrium shape, specific proper- 
ties for excitation energies and electromagnetic transition probabilities show up. Due to this 
reason, one customarily says that static values of intrinsic coordinates determine a phase 
for the nuclear system. A weak point of the boson description with a complex anharmonic 
Hamiltonian consists of the large number of the structure parameters which are to be fitted. 
A much smaller number of parameters is used by the coherent state model (GSM) [5] which 
uses a restricted collective space generated through angular momentum projection by three 
deformed orthogonal functions of coherent type. The model is able to describe in a realistic 
fashion transitional and well deformed nuclei of various shapes including states of high and 
very high angular momentum. Various extensions to include other degrees of freedom like 
isospin ^, single particle [3] or octupole 8| degrees of freedom have been formulated ol. 

It has been noticed that a given nuclear phase may be associated to a certain symmetry. 
Hence, its properties may be described with the help of the irreducible representation of the 
respective symmetry group. Thus, the gamma unstable nuclei can be described by the 0(6) 



symmetry 



lOj, the symmetric rotor by the SU{3) symmetry and the spherical vibrator 



by the U{5) symmetry. The gamma triaxial nuclei are characterized by the invariance 



symmetry group D2 



ll| of the rigid triaxial rotor Hamiltonian. Thus, even in the 50's, 



the symmetry properties have been greatly appreciated. However, a big push forward was 
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brought by the interacting boson approximation (IBA) , Q , which succeeded to describe 
the basic properties of a large number of nuclei in terms of the symmetries associated to 
a system of quadrupole (d) and monopole (s) bosons which generate a U{6) algebra. The 
three limiting symmetries U{5), 0(6), SU{3), are dynamic symmetries for U{6). Moreover, 
for each of these symmetries a specific group reduction chain provides the quantum numbers 
characterizing the states, which are suitable for a certain region of nuclei. Besides the virtue 
of unifying the group theoretical descriptions of nuclei exhibiting different symmetries, the 
procedure defines very simple reference pictures for the limiting cases. For nuclei lying close 
to the region characterized by a certain symmetry, the perturbative corrections are to be 
included. 

In Ref. [l^, a new classification scheme was provided, all nuclei being distributed within a 
symmetry triangle. The vertices of this triangle symbolize the U{5) (vibrator), 0(6) (gamma 
soft) and SU{3) (symmetric rotor), while the legs of the triangle denote the transitional 
region. Properties of nuclei lying far from vertices are difficult to be explained since the 
states have some characteristics of one vertex while some others are easy to be described by 
using the adjacent symmetry. The transition from one phase to another reaches a critical 



point depending on the specific parametrization as well as on the transition type 



1. In Ref. 



18 



16 



17? 



19| . it has been proved that on the U{5) — 0(6) transition leg there exists 



a critical point for a second order phase transition while the U{5) — SU{3) leg has a first 
order phase transition. 



Recently, lachello [20|, l2l| pointed out that these critical points correspond to distinct 
symmetries, namely -^(5) and X{5), respectively. For the critical value of an ordering 
parameter, energies are given by the zeros of a Bessel function of half integer and irrational 
indices, respectively. The description of low lying states in terms of Bessel functions was 
used first by Jean and Willet [10|, but the interesting feature saying that this is a critical 
picture in a phase transition and defines a new symmetry, was indeed advanced first in 
Ref.fl. 

Representatives for the two symmetries have been experimentally identified. To give 



an example, the relevant data for ^'^^Ba [22] and ^^^Sm 23| suggest that they are close 



to the -£^(5) and X{5) symmetries, respectively. Another candidate for E(5) symmetry. 



26|, the 



proposed by Zamfir et al. [25] is ^°^Pd. Using a simple IBA Hamiltonian, in Ref. 
low lying spectrum of ^^'^Pd is realistically described. Comparing the E(5) predictions with 
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the experimental data concerning energy ratios in the ground band and the normahzed E2 
transition probabihties for the states 4"*" and 0^, one concludes that this nucleus is a good 
E(5) candidate. However, in order to decide which Pd isotope is closer to an E(5) behavior, 
further investigations are necessary. A systematic search for E{5) behavior in nuclei has 
been reported in Ref. 24l |. 

Short after the pioneering papers concerning critical point symmetries appeared, some 



other attempts have been performed, using other potentials like Coulomb, Kratzer [27| and 
Davidson potentials [28|. These potentials yield also Schrodinger solvable equations and the 
correspon ding results may be interpreted in terms of symmetry groups. 



In Ref. |30l| we advanced the hypothesis that the critical point in a phase transition is 
state dependent. We tested this with a hybrid model for ^^^Ba and ^°'^Ru. 

The departure from the gamma unstable picture has been treated by several authors 
whose contributions are reviewed by Fortunato in Ref. 33|. The difficulty in treating the 
gamma degree of freedom consists in the fact that this is coupled to the rotational degrees of 



freedom. A full solution for the Bohr-Mottelson 
of gamma deformation can be found in Refs. 34 



iamiltonian including an explicit treatment 



35| . Therein, we treated separately also the 



gamma unstable and the rotor Hamiltonian. A more complete study of the rotor Hamiltonian 
and the distinct phases associated to a tilted moving rotor is given in Ref. 36] . Distinct 
solutions, expressed in laboratory frame shape coordinates, have been reported in Refs. 



|37|, 



38 



4J]. The gamma dependent part of the wave function has been found as a solution 



of a specific differential equation in Ref. 4ol |. 

Finding the 7 depending part of the wave function becomes even more complicated when 
we add to the liquid drop Hamiltonian a potential depending on (3 and 7 at a time. To 
simplify the starting problem related to the inclusion of 7 one uses model potentials which 
are sums of a /3 term, V{P), and a factorized /3 — 7 term f/(7)//5^. In this way the nice 
feature for the beta variable to be decoupled from the remaining four variables, specific to 
the harmonic liquid drop, is preserved. In the next step, the potential in gamma is expanded 
around 7 = or 7 = |. In the first case if only the singular term is retained one obtains 
the infinite square well model described by Bessel functions in 7. If the 7^ term is added 
to this term, the Laguerre functions are the eigenstates of the approximated 7 depending 
Hamiltonian, which results in defining the so called X(5) approach. 

Note that any approximation applied to the 7-Hamiltonian modifies automatically the 
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differential equation for (3. Indeed, the centrifugal term r(r + 3)//3^ disappears but another 
one is expected to come from the /3 — 7 coupling after some approximations are performed. 

The drawback of these approximations consist in that the resulting functions are not 
periodic, as the starting Hamiltonian is. Moreover, they are orthonormalized on unbound 
intervals although the underlying equation was derived under the condition of I7I small. 
Moreover, the scalar product of the resulting functions is not defined with the integration 
measure |sin37|(i7 as happens in the liquid drop model. Under these circumstances it 
happens that the approximated Hamiltonian in 7 looses its hermiticity. 

In a previous short publication 4l|] we proposed a scheme where the 7 Hamiltonian 
is solvable Hamiltonian. Moreover, its eigenstates are the spheroidal functions, which are 
periodic. Also one proves that the X(5) model is recovered in the limit of 7-small. Here we 
give details about the calculations and describe some new numerical applications. Moreover, 
we complete our formalism by considering a Schrodinger equation in the beta coordinate, 
involving a potential of the Davidson's type. 

By the numerical applications we want to see whether curing the mentioned drawbacks 
for gamma wave functions and Hamiltonian, would bring substantial quantitative corrections 
to energies and E2 transitions given by the X(5) formalism. Even if the found corrections 
are not dramatic the virtue of the proposed formalism to extend the X(5) description to 
a theory which is consistent with the symmetries of the starting Hamiltonian as well with 
some basic principles of Quantum Mechanics concerning the hermiticity property of the 
model Hamiltonian remains an important achievement of the present paper. 

The above objectives are reached according to the following plan. The starting Hamil- 
tonian is presented in Section II, where the separability conditions are discussed. Several 
methods for treating (3 are described in Subsection A while the Hamiltonian in 7 is studied 
in Subsection B. The specific procedure of the present paper is presented in Section HI. 
Numerical applications are given in Section IV. The final conclusions are summarized in 
Section V. 
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II. THE STARTING HAMILTONIAN 



Written in the intrinsic frame of reference, the original Bohr-Mottelson Hamiltonian has 
the expression: 



2B 



I d 1 d . ^ d__J_ ^ Ql 

(3^ d(3^ d(3 /?2 sin 37 d-f ^'"^ ^d-f A(3^ ^^^,3 sin'(7 - l^^k) 



(2.1) 



where the dynamic deformation variables are denoted by (3 and 7 while the intrinsic angular 
momentum components by Qk-, with k = 1,2, 3. Within the liquid drop model the potential 
energy depends quadratically on /5. Here we assume that the potential energy depends 
on both deformation variables, beta and gamma. Without exception the solvable models 
proposed for a simultaneous description of /? and 7 variables adopt the variable separation 
methods. Two situations are to be distinguished: 

a) If the potential energy term is depending on deformation variables in a separable 
manner: 

V{P,^) = V{P) + U{^), (2.2) 

and some additional assumptions are adopted, the eigenvalue equation associated to H (12. ip 
can be separated in two parts, one equation describing the beta variable and the other one 
the gamma deformation and the Euler angles Q = {61,62,63). Indeed, in this case the 



separation of variables, achieved by various models, is based on two approximations j58|, |59| : 
i) restriction to small values of 7, i. e. I7I ^ 1; ii) replacing the factor 1//?^ by 
in the terms involved in the equation for 7. The diagonalization of the Bohr-Mottelson 
Hamiltonian shows that the first approximation is valid for large 7 stiffness while the second 



one for small 7 stiffness 



59] 



b) A complete separation of equations for the two variables, (3 and 7 is possible if we 
choose the potential 

\/(/3,7) = m + f/(7)//9'. (2.3) 

In this way the approximation consisting of replacing jS"^ by {(3^) is avoided but the restriction 
to the case I7I <^ 1 is still kept. Thus, the theory involves one parameter which is the 7 
stiffness which affects the excitation energies in both the beta and gamma bands. 
In what follows the two equations, for P and 7, will be considered separately. 
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A. The treatment of the /? Hamiltonian 



The solvable models for /3, presented here, have been used by E(5) formalisms, which 
ignore the potential in 7. Considering the potential in 7, of course, the picture for (3 is 
changed. However, as we shall see later on, it is very easy to derive analytically the energies 
and wave functions associated to (3 from the corresponding results of the E(5) descriptions. 
Actually, this is the motivation for reviewing, here, the beta solvable models. 

The equation in /? is: 



1 a ^,8 



+ 



A 



+ u{(3) 



(2.4) 



(3^8(3"^ d(3 /32 

where A is the eigenvalue of the Casimir operator of the S0{5) group. This is related with 
the seniority quantum number r, by A = r(r + 3). The 'reduced' potential u{P) and energy 
e are defined as: 

E^^e, V^-u. (2.5) 

where E denotes the eigenvalue of the Hamiltonian H corresponding to the potential V{13). 
Here we mention the most used potentials for (3: 



1. The case of u{p) = 



A full description of the eigenstates of the Bohr-Mottelson Harniltonian satisfying the 
symmetry f/(5) D SO (5) D SO (3) D 50(2), may be found in Refs.[34]. In particular, the 
solution of the radial equation (12.41) with u{P) = is easily obtained by bringing first Eq.( 
12. 4p to the standard Schrodinger form by changing the function f to by: 



The equation obeyed by the new function ip, is: 



+ 



{T + l){r + 2) 
P' 



^ = 0. 



This equation is analytically solvable. The solution is: 

i^nriP) = , T.r-^^^^^^^J;'''/'(/5')/5^"''exp(-/5V2), 



\ r(n + r + 5/2)' 
en = 2n + T + 5/2, n = 0, 1, 2, r = 0, 1, 2, 3, 



(2.6) 



(2.7) 



(2.8) 
(2.9) 
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where LJ^ denotes the generahzed Laguerre polynomials. The number of polynomial nodes 
is denoted by n and is related to the number of the quadrupole bosons (A^) in the state, by: 
N = 2n + T. Consequently, the initial equation fl2.4p has the solution 



friT P '4^nT ■ 



(2.10) 



The spectrum, given by Eq. (12. 9p . may be also obtained by using the unitary represen- 
tation of the SU{1, 1) group with the Bargman index k = (r + 5/2)/2. Indeed, the standard 
generator for SU{1, 1) are: 



where 



^^,^(r+l)(r + 2) 



13' 



Hq obeys the following equations: 



(2.11) 



(2.12) 



(2.13) 



2. Davidson's potential 



Another potential in f3 which yields a solvable model is due to Davidson [42|]: 



2 I Po_ 



(2.14) 



This potential has been used by several authors in different contexts 28|, |43|, |4J]. The 
potential has been used by Bonatsos at al 59|], to describe the dynamic deformation variable 
p. For this potential the above equations (2.6-2.10) hold for r replaced [28i] by 

1/2 



(2.15) 



In particular, the excitation energies have the expressions: 



2n + l 



3\2 



1/2 



(2.16) 
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The factor (3^ is considered to be a free parameter which is to be determined variationally 
for each angular momentum, as suggested in Ref. 



df3i 



0, 



(2.17) 



where R denotes the ratio of the excitation energy of the ground band state L"*" and the 
excitation energy of the state 2+. 



3. Five dimensional infinite well 

Now, let us turn our attention to the situation considered by lachello in Ref.f^ol, where 
the potential term associated to the spherical to gamma unstable shape transition is so flat 
that it can be mocked up as a infinity square well 

'0, P< (3^, 

, OO, (3 > Py,. 

A more convenient form for the equation in /3, is obtained through the function transforma- 
tion: 

m = p"'fm. (2.19) 



u{(3) 



(2.18) 



The equation for ip is 



(pLf 1 dip 
d(P^]3d(3^ 



f = 0. 



Changing the variable (3 to z hj 

z = kp, k = i/e 

and denoting with (p{z) = if {(3) the function of the new variable, one arrives at: 

[r + 3/2)2- 



d'^if 1 dip 
dz"^ z dz 



1 



ip = 0. 



(2.20) 



(2.21) 



(2.22) 



This equation is analytically solvable, the solutions being the Bessel functions of half integer 
order, Jt+3/2{z). Since for [3 > j3yj the function (p is equal to zero, the continuity condition 
requires that the solution inside the well must vanish for the value of (3 equal to (3yj. This, 
in fact, yields a quantized form for the eigenvalue E. Indeed, let x^^r be the zeros of the 



Bessel function J,, : 



^r+3/2(a:^,r) = 0, ^ = 1, 2, . . . ; T = 0, 1 , 2, 



(2.23) 



Then, due to the substitution introduced in Eqs.(2.21) and (2.5) one obtains: 

E,. - ^*|.. - I;. (2.24) 

Concluding, the differential equation for the beta deformation corresponding to an infinite 
well potential provides the energy spectrum given by Eq. fl2.24p and the wave functions: 

/e,. = C5,./?-3/V.+3/2(f^/3), (2.25) 

where C^^r is a normalization factor. 

It is worth noticing that the spectra corresponding to E(5) and Davidson potentials 
become directly comparable by establishing the formal correspondence n = ^ — 1. 



4-. A hybrid model 

In ref . jsO] we advanced the idea that the critical point for a phase transition is depending 
on the nuclear state. Therefore the system may reach the critical point in a state of angular 
momentum J, but in a less excited state, like (J — 2) + , the system could behave according 
to the initial nuclear phase. 



According to Ref. 



301 ] the potential energy in the beta variable is depending on angular 



momentum in the following way: 



r if < /? < oo, L < 2, 
u{f3) = I 0, if < /3 < L > 4, 

^ oo, if < /5 < oo, L > 4. 

The states of interest and their energies have the following expressions: 



\LtM) 



En 



2n\ 



r(n + r + 5/2) 



/?^Lr^/^(/5^)e-/^^/^G^f(7,fi), 



2B 



(2n + r + 5/2), (n, r) = (0, 0), (0, 1), L = 2r, 



|L+,M) = Q,./3-3/V.+3/2(/3x^,.//5^„)Gf^{,.(7,f^), 



E, 



2B PI 



'2 ' 
w 



(e,r) = (l,2), (1,3), (2,0). 



(2.26) 



(2.27) 



The factor functions depending on the beta variable are solutions of Eq. fl2.4p with the 
reduced potential given by Eq. (l2.26p . The equation for 7 deformation and Eulerian angles 
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(O) has the solution G^^ . A possible excited state phase transition was pointed out in 



Ref. 



3l| . by using a sixth order solvable boson Hamiltonian. A potential of an intrinsic 



deformation radial variable r, involving a centrifugal term and a + term, shows up. In 
the state lying at the top of the barrier separating the two wells of the potential, the system 
undergoes a phase transition. This issue has been also addressed by studying the energies 
and wave functions singularities for Lipkin model as well as for a two level boson model 32|. 

Note that in all treatments mentioned above, no potential in 7 is considered. Due to this 
fact the spectra and wave functions are labeled by the seniority quantum number r. This 
feature does not hold when we switch on the 7-depending potential and moreover impose 
variable separability by approximating the terms depending on 7 



B. The description of 7 degree of freedom 

Let us consider the Hamiltonian 

H = -—^jr- sin 37^ + U{^) + W{^, Q), (2.28) 
where f/ is a periodic function in 7 with the period equal to 27r and 

WbM = \t-^^^Ql (2.29) 

with Qk denoting the components of the intrinsic angular momentum. 



1. Violating some basic properties 

Any approximation for the potential, by expanding it in power series of 7, alters the 
periodic behavior of the eigenfunction. Moreover, the approximating Hamiltonian loses its 
hermiticity with respect to the scalar product defined with the measure for the gamma 
variable, | sin(37)|(i7. 

We illustrate this by considering the case of a little more complex potential 

= Micos(37) + U2Cos2(37). (2.30) 

Performing the change of function ip = y^|sin(37)|'?/', the eigenvalue equation Hip = Eip, 
becomes Hip = 0, with 
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H 



&^ 9 
+ 4 



1 + 



1 



sin (37) 



-U -W + E. 



(2.31) 



We shall consider two situations: 

A. Suppose that I7I -C 1. Expanding the terms in 7 in power series up to the fourth 



order, one obtains: 



W4 



Ml + U2 
1 



97'(y + «2 



277' 



Ml 



U2 



-11 + 27^ + ^ 

3 I ^ 9 



(2.32) 



4 I 72 ^ 3 15 189 



Ql 



The low index of U and W suggests that the expansions in 7 were truncated at the fourth 
order. Details about the approximations involved in the following derivation may be found 
in Appendix A. 

Note that due to the term W, the equations of motion for the variable 7 and and Euler 
angles are coupled together. Such a coupling term can in principle be handled as we did for 
the harmonic liquid drop in ref. 34, |35l]. Here, we separate the equation for 7 by averaging 

with an eigenfunction for the intrinsic angular momentum squared. The final result for 

is: 



= ^ + ^(l-(Ql))+/^o + W + /^47^ 

+^ (1 + 27^) {Ql - Ql) , (2.33) 

ho = E- ^L{L + 1) + ^ (Ql) - {Uo + Ml + U2) + y , 

2 13 / \ 9 27 

h2 = --L{L + 1) - ^ {Ql) + + 9m2 + -, 

/14 = --L{L + 1) - — (Ql) - -m - 27m2 + -. 
27 ^ ' 126 \^^/ 8 14 

where L denotes the angular momentum. If the average is made with the Wigner function 
D^j^, important simplifications are obtained since the following relations hold: 

{Ql-Ql) = 0, {Ql) = K' (2.34) 
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Actually, this is the situation considered in the present paper. Note that H4 contains a 
singular term in 7, at 7 = 0, coming from the term coupling the intrinsic variable 7 with 
the Euler angles. One could get rid of such a coupling term by starting with a potential in 
gamma containing a singular term which cancels the contribution produced by the W term. 
Thus, the new potential would be 

U' = U+ \, (2.35) 
4sin2(37) ^ ' 

The corresponding fourth order expansion for the Hamiltonian is: 



= ^ + ^ + K + Kl' + Kl\ (2.36) 

27 27 

K = ho + K\ h'^ = h2 + —K^, K = K + —K\ 

Some remarks concerning the equation H'^ip = are worth to be mentioned: 

i) If in this equation one ignores the 7^ term, the resulting equation has the Laguerre 
functions as solutions and moreover the Hamiltonian exhibits the X{5) features. 



ii) Note also that the Hamiltonian coefficients are different from those of Ref. 
difference is caused by the fact that here, the expansion is complete. 



33|. The 



iii) Taking in the expanded potential ui = U2 = and ignoring, for 7 small, the term 

20-' 



^^r/r^7^, the resulting potential is that of an infinite square well which was treated by lachello 



in Ref 2l| . The solutions are, of course, the Bessel functions of half integer indices. 

iv) Irrespective of the potential in 7, in the regime of I7I small a term proportional to 7^ 
shows up due to the rotational Hamiltonian W. Therefore, even in the case the potential 
is taken as an infinite square well, of the form 1/7^, the equation describing the 7 variable 
admits a Laguerre function as solution and not, as might be expected, a Bessel function of 
semi-integer index. Amazingly, the potential in 7 is also of Davidson type. 

v) None of the mentioned solutions is periodic. 

vi) Also the approximated Hamiltonians are not Hermitian in the Hilbert space of func- 
tions in gamma with the integration measure as introduced by the liquid drop model, i.e. 

I sin S'yld'j. 

B. The case I7 — 7r/6| ^ 1. Using the fourth order expansion in y = I7 — 7r/6|, given in 
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Appendix A, one obtains a Hamiltonian similar to that given by Eq. (2.36): 

H'. = ^,+ Kl' + + 2V?>y [l + ^y^) {Ql - Ql). (2.37) 

If {Q\ — Q2) = and, moreover, one ignores the term in 7^ the resulting equation in 7 
describes a harmonic oscillator. Again the eigenf unctions, i.e. the Hermite functions, are 
orthogonal on an unbound interval of 7, and not on [0, 2tt]. 



2. Toward an exact treatment which preserves periodicity and hermiticity 

In order to remove the drawbacks mentioned above, we try first to avoid making approx- 
imations. Thus, let us consider the Hamiltonian given by Eq. f l2.28p where instead of U we 
consider U' as defined by Eq. (12.351) . and ignore for a moment W . Changing the variable 
X = cos 37, the eigenvalue equation associated to this Hamiltonian, HS = ES, becomes: 

Note that we denoted the eigenfunction by S which suggests that the differential equation 
(I2.38P is obeyed by a spheroidal function. If ui = U2 = K = 0, the solution of this equation is 



the Legendre polynomial P„ while E = 9n{n + l). This case has been considered in Ref. |33|. 
This function may be used to approximate the solution of the original liquid drop model. 
For other particular choices of the coefficients ui,U2 defining the potential in gamma, the 
solution is readily obtained if one compares the above equation with that characterizing the 
spheroidal oblate functions [45|] 

(1 _ ^ _ + f A_ - cV - S^m = 0. (2.39) 

^ ^ dx^ ax \ 1— x^y 

The prolate case is reached by changing c ic. 

For c = 0, the solutions of Eq. (12.391) are the associated Legendre functions P™. For c 7^ 0, 
Snmi with m, n integers and n > m > 0, are linear series of these functions. 

In the case Mi = 0, the solution of Eq. (l2.38p is identified as being the spheroidal function 
while the energy is simply related to Xnm'- 

2 ^2 , 1 

m = C = — , Xnm = -Enm- (2.40) 

z y y 

Here E^m denotes the eigenvalue E corresponding to the quantum numbers n and m. 
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FIG. 1: (Color online). The spheroidal energy E'^^^ = A„m = Enm/^, for < m < n < 3 are plotted 
as functions of c = y^/S. 

For |c| small the energies Enm exhibits the asymptotic expansion 

2 (n(n + 1) + - 1) 

Enm ~ 9n n + 1 V tt r^^s 2.41 

^ ^ (2n-l)(2n + 3) ^ ' 

1 [{n — 1)^ — m^] (n^ — m^) 2 

^ 18 (2n - 3)(2r2 - l)3(2n + 1)^^ 

1 [{n + 1)2 - m2] [{n + 2)^ - m^] ^ 



18 (2n+ l)(2n + 3)3(2n + 5) 
Eg. (12.411) considered for a fixed m but various n, defines a band. Similar expansions may 
be derived for |c| large. 



Enm ~ -U2 + 3gv^ + 9 ( - 1^ I - -^%(11 + - 32m2), 

V 8 / QAJu^ 



q = 2(n - m) + 1 (2.42) 

We remark that the spectrum has a rotational behavior for small c, due to the term n{n + l) 
while for large values of c it has an oscillator feature, the energy depending linearly on n. 
If one needs the expansion up to the l/c^ terms, the results for the first few energies are: 

1 2 5 33 \ 
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9 69 417 \ , , 

It is worth spending few words about Fig. 1 where the energies correspond to the 
spheroidal functions with the parameters specified by Eq.(2.40). Indeed for c ^ one 
notices some multiplet degeneracy which suggest a symmetry with respect to K, i.e. a ro- 
tation invariance of states of a given n. Increasing c the split in energy is similar to that in 



Nilsson 



60| model when the energy is Q dependent. The difference is that while in Nilsson 



model each deformed state is a superposition of states with different angular momentum, 
here the multiplet members are characterized by the same n. In this respect the feature 
shown in Fig.l is similar to the one obtained with a spherical projected single particle basis 
6l| . In the region of large c, for a given large n the set of states of different m seem to 
form a band. On the other hand for a fixed m the set of states with different n is a band of 
equidistant energy levels. 



C. Approximation which does not affect periodicity and hermiticity 

Now, we shall focus on an approximate solution which preserves the periodicity in 7. For 
that purpose we consider the Hamiltonian 

H = ———sm3^— + U{^), 
sm 37 c'7 0*7 

UM = Ml cos 37 + M2 cos^ 37 H 5—. (2.44) 

4 sin 7 

Changing the function by the transformation \E' = | sin(37)|^^/^$, for sin(37) 7^ 0, the 
eigenvalue equation for H is if$ = with H given by Eq. (2.31) for W = 0. 

Under the regime of I7I small, we take the 0(7^) expansion of the terms depending on 
7 and in the final expression approximate 7 ~ sin 7. In this way the eigenvalue equation 
becomes: 
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1 n „ + a - 2g COS 27 5— $ = 0, with 

yc/7"' 4 sin 7 y 

1 9 9 347 ^ 10 

q = - + -Ui + -M2, u = U2 + -7-7, a = + — g + u. 



(2.45) 



3 8 4' 108' 9 

We suppose now that this equation is vahd in the interval [0,27r]. The equation f l2.45p is 
just the trigonometric form of the spheroidal functions. The algebraic version is obtained 
by changing the variable x = cos 7. 

For K = 1 one obtains the Mathieu equation: 



+ a - 2gcos27 $ = 0. 



(2.46) 



There are two sets of solutions, one even and one odd denoted by $^(a, q, 7) and $ (a, q, 7), 
respectively. For q = 0, both solutions are periodic for any positive value of a. 



$+(a, 0, 7) = cos [y/aj'j , $ (a, 0, 7) = sin (v^7) . 



[2.47) 




10 15 



FIG. 2: (Color online)The characteristic curves are plotted as functions of q for several values 
of n. 

For q ^ the Mathieu functions are periodic in 7 only for a certain set of values of a, 
called characteristic values. These are denoted by for even and c~ for odd functions, 
respectively. In the plane (a,g), the characteristics curves separate the stability regions, 
shown in Fig. 2 by gray color, from the non-stability ones, indicated by white color in 
the quoted figure. For g = the equalities c^(0) = ri^ hold. By means of Eq. (12.451) the 
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characteristic values determine the energy E. Thus, the energy spectrum is given by 



u 



with Ef^ = — yg. The corresponding wave functions are the eUiptic cosine and eUiptic 
sine functions respectively: 

*o(?>7) = ^=ceo(g,7),*n(?,7) = ^ce„(g,7), 
V27r 



7) = ^se„(g, 7), n = 1, 2, ... (2.48) 



They form an orthogonal set. The matrix elements of the gamma depending factors of the 
transition operator can be easily calculated in Mathematica. Moreover, in the regime of |g|- 
small these matrix elements can be analytically performed, since the following representation 
of the wave functions hold: 



$^(7) ^ cos(?27 - 9±) 



cos [(n + 2)7 - e±] cos [{n - 2)7 - e±] 



q\ (2.49) 



4(n + l) 4(n-l) 
where n > 3, ^+ = and = n j1. The corresponding energies have very simple expressions: 



E^ ^ u — —q 



2 

9 " 2'' 



Et ^ E^ ^u - —q - J, 

10 

E^ ^ E^ ^u + A- —q-^.. 



10 0^ 

^ + J_ n>3. (2.50) 



Normalizing the above functions to unity, on the 7 interval [— vr, tt] with respect to the 
integration measure ^7 and calculating with the resulting functions the matrix elements of 
the 7 depending factors involved in the electric transition operator one obtains the curves 
represented in Fig. 3. 

Obviously, a phase transition is determined by the combined effects coming from the 
behavior of the wave function in the f3 and 7 variables, respectively. 

In the X{5) formalism, the eigenfunction of the (3 Hamiltonian is a Bessel function of 
irrational index, while the 7 Hamiltonian's eigenfunction is a Laguerre polynomial. 

Here we propose to change the description in the 7 space either by a spheroidal or 
by a Mathieu function. These functions are periodic and the corresponding Hamiltonians 
Hermitian. Moreover, in both versions, the X{5) Hamiltonian is obtained in the limit of 



small I7I. In Ref. 46] a periodic 7 potential with a minimum in 70 = 7r/6, was considered. 
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FIG. 3: The matrix elements of T^^, for cos 7 (ee' = +) and sin 7 (ee' = — ) are represented as 
functions of q. 

The model is solvable and the wave function is a Legendre polynomial. Moreover, energies 
are analytically obtained. By contrast the situation considered here is more complex, the 
energies being obtained by solving the eigenvalue equation for the spheroidal functions. 



However the picture described in Ref. 



46| is recovered under some particular restrictions. 



D. Including the rotational term preserves periodicity and hermiticity 

We recall that so far the rotational term W was left out. Now we turn our attention 
to this term. If we average W with the Wigner function -Df/^ and add the result to the 
potential U' given by Eq.(2.35) and then following the same path as before, one ends up 
also with an equation for a spheroidal function. Indeed, let us consider the average of W: 

{LK\W\LK) = - , ^\ 

8 sm 37 8 sm 7 

D = L{L + l)-K^-2. (2.51) 

When I7I <^ 1, this expression admits the following second order expansion in sin 7: 

^(7) = -r-^ + o + 1) - - 2! (1 + 2 sin^ 7) . (2.52) 
4 sm 731- J ^ ^ 

The term L{L + l)/3 from the above expression, multiplied with the factor 1//?^ are added 
to the equation describing the variable (3. In the case one makes the option for a infinite 
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well potential in /3, the renormalization just mentioned leads to an equation in /3, whose 
solution is a Bessel function with the index 

z.= + + (2.53) 

For the case I7I <^ 1, we consider the second order expansion in sin 7 for the full Hamil- 
tonian (2.28). The result is a trigonometric form for the equation of the spheroidal function: 

§^+(e'- - C sin^ 7) = 0, (2.54) 
0*7^ Y 4 sm 7 J 

with the notations: 

q D 2D 

E' = E + --UI-U2--, C = — - — -9U2. (2.55) 

4 3 3 2 ^ ^ 

As suggested by the expression of the starting Hamiltonian, the remaining terms of ^^(7) 
should be multiplied with This coupling of (3 and 7 variables is usually considered as 

a renormalization term for the potential in 7, by replacing the factor 1//?^ by the constant 
1/ The notation is used for the expectation value of in the ground state which 
results in redefining the constant C, in Eq. (2.55). 

Eq. (12.541) can be brought to the form given by Eq.(2.39) by making a successive change 
of function S = \sm'y\)~^^'^ cp and variable, x = cos 7. Indeed, the resulting equation is that 
of the spheroidal function defined by: 

7 1 

Xnm = Enm + -U^ + 8U2 + 2 - D + -L{L + 1), 

= \ui + 9m2 - \d. (2.56) 

This equation has been used by some authors of this paper, in Ref. |4l|. to describe the 
spectrum and the E2 properties of ^^^Sm. 



E. Recovering X(5) in the limit of |7|-small 



It is worth comparing the present formalism based on the spheroidal functions with the 
X(5) approach. It is easy to prove that, indeed, X(5) is the limiting case for our approach. 
Indeed, considering the second order expansion in 7 of the terms involved in Eq. (12.541) one 
arrives at: 



E' 



1 



4 



1 



60 



0. 



(2.57) 
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This equation is characterizing the X{5) model, with all harmonic contributions included. 
Indeed, changing the variable ^ = 57, with q — \JC -\- the differential equation be- 



comes: 



+ 



1 /^/ 1 



l\ 1 



47 e 



(f, a 



K 



(2.58) 



Comparing this equation with that describing an harmonic, isotropic plane oscillator: 



de 



+ 



2(n + a + l) -C' 



a 



4) e 

one identifies the function (p with the function Laguerre: 



^ r 



K 

0, a = —. 



while the system energy is: 



— ,cLi{e) exp (--e 

r(n + Q; + l)^ ^ ^ V 2^ 



Er, 



2 C + 



K , 

n+ — + l +'Ui + 'U2 + -tK' 



3. 



(2.59) 



(2.60) 



(2.61) 



60 y V 2 J 4 

The property of reaching the X{5) model in the hmit of small values of I7I, holds also for 
the Mathieu functions. Indeed, these functions satisfy a differential equation which is of 
spheroidal type. Consequently, in the limit I7I ^ 1, the Mathieu functions may account 
for the properties which are specific to the -^(5) approach. Numerical applications with 
Mathieu functions will be published elsewhere. 



III. THE PRESENT APPROACH 



Here we summarize the procedure adopted in the present paper, to treat a phenomeno- 
logical solvable Hamiltonian defined in the space of the variables f3 and 7. 

The potential in the two variables is considered to be of the form given by Eq. (2.3). As 
V{(3) we take the Davidson potential (2.14). Including the terms proportional to (this 
is 3^ [L{L + 1)]) from the rotational term in the Schrodinger equation for /3 , the resulting 
equation admits solutions which formally coincide with those given by Eqs. (2.8) and (2.9) 
but having instead of r a irrational quantum number p defined as: 



1/2 



(3.1) 



This expression is obtained by writing the coefficient of from the Schrodinger equation 

associated to the (3 variable in the form: 



(p + l)(p + 2) = 2 + -L(L + l)+/3ol 



(3.2) 
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This equation has two solutions, one written above, while the second one is differing from 
the first one by the sign of the square root term. Let us denote for a while the two solutions 
by pj-. Note that the Davidson potential is not a continuous function in /3 = 0. This causes 
the fact that for L=0, we have: 

lim p. = 0, lim p_ = —3. (3.3) 
Consequently, the corresponding spectra are given by: 

Ei^^ =2^ + ^' 

Ei^-) = 2n-3 + ^. (3.4) 
Therefore, the full spectrum of the 5-dimensional oscillator is recovered only if both solutions 



are considered at a time 62] • Since, as we shall see a bit later, /Sq is far from origin we make 
the option for the branch corresponding to the eigenvalue p+. The choice is justified by the 
fact that for /3o 7^ the wave function corresponding to p- is singular in origin and therefore 
it is not a convenient solution. 

As explained before, for the states belonging to the ground band Pq was fixed variationally, 
by Eq.(2.17). We extended Eq. (2.17) to the beta and gamma bands, respectively. The first 
derivatives of the ratios R^i \ defined by 

(3) - Eo+ 

RT = ^ 7^.^> 4, 

L, 7-1 Z? ' — ' 

-C/9+ — -C/n+ 

= '\ L> 4, (3.5) 

have the (3q dependence shown in Fig.4 for some particular values of L. We fix (3q for the 
states in the band k {=(3, 7) so that the first derivatives of the ratios R^^'' are maximum. 
From Fig. 4 one sees that each curve has a well pronounced maximum. Collecting the values 
of obtained in this way, and representing them as function of L one obtains a straight 
line for both beta and gamma bands, as shown in Fig. 5. Extrapolating these straight lines 
for L = 0, 2 in beta band and L = 2, 3 in the gamma band, one obtains a one to one 
correspondence between the states in the two bands and the values of Pq. 
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FIG. 4: The first derivative of the ratios Rj^' FIG. 5: The solutions ' of the equation 
(upper panel) and R^?^ (lower panel), defined by ^^^2— = for A; = /3 (upper panel) and A; = 7 
Eq. (2.5), are plotted as function of /Jq. (lower panel), arc plotted as functions of the an- 

gular momentum L. 
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Let us denote by 2n + p + |) the energy provided by the Shrodinger equation 

associated to the variable p. The index D suggests that the P potential is chosen to be of 
Davidson's type. 

For comparison we considered also an infinite square well potential for the j3 variable. In 
this case the energies associated to the (3 variable are denoted by E^g''^\ They correspond 
to the Bessel function of index s + 3/2 with s defined by the following equation: 

(5+2) =4 + 3^(^ + 1) (3.6) 

^ is an ordering index for the zeros of the Bessel function. Therefore the irrational index for 
the Bessel function will be: 



« = + + (3.7) 

Other possible ways of renormalizing the differential equation for /3 are discussed in Appendix 

C. The cases when the spheroidal function formalism for the 7 variable is considered at a 
time with an oscillator potential or a hybrid potential potential in (3 will discussed elsewhere. 
As for the potential in 7 we considered 

1 



f/(7) 



9 



Ui cos 37 + U2 COS 37 + 



(3.8) 



4 sin^ 37 

Assuming that I7I ^ 1 the rotational term is expanded in powers of sin 37. From the 
rotational term we depict the term not depending on 7 and proportional to L{L + 1), 
otherwise being proportional to and add it to the Hamiltonian in /? which results in 
having a renormalization of the centrifugal term. In the remaining terms, we approximate 
by Thus, the Hamiltonian in 7 will comprise an overall factor If in the 

Hamiltonian which multiplies this factor, one changes the function ip ^ S — \ sin37|~^/^<^ 
and the variable 7 — > a; = cos 37, the corresponding Schrodinger equation is that of a 
spheroidal function defined by Eq. (2.39) with 



A„„ = i _ 1 -^D + h{L + 1)) , 
2 1 /I 2 \ 

C = - -Ml + Uo D , 

9 V2 ' 27 ; ' 

m = -. (3.9) 
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For illustration, in Figs. 6, 7 we give few potentials corresponding to {j3) = 1 and different 
sets of (mi,M2). The spheroidal functions normalized to unity in the interval [0,7r/3], given 
by Eq. (2.39) with the parameter determined by U2) used in Figs. 6, 7 are represented 
as functions of 7 in Figs. 8, 9 for three L levels from the ground band. Once we fix the 7 
potential, we can calculate the energy associated to the 7 variable. 

The total energy for the system described by the decoupled P and 7 variables is: 

Et]n'm;LK = Eo + AE^l^'^ + FEljl k = D,B. (3.10) 

The coefficient F includes the factor 1/ mentioned above. Due to this feature there is no 
need to specify the average value of Note that the energy determined by the rotational 
degrees of freedom has been already included when the term of the Hamiltonian denoted by 
W was averaged with the Wigner functions. 

If in the space of (3, the Schrodinger equation involves the Davidson's potential, the wave 
function describing the whole system is: 

\np- n'm- LMK) = m^^{(3)Sn'm{cos^l)^^^^^ {d^k + (-l)''+'^/^if,_i^) ,^ = y- 

(3.11) 

The ground, beta and gamma bands are defined by the quantum numbers: 

n = 0, n' = 0, m = 0, K = 0, L = 0,2, ... ground band, 
n = 0, n' = 1, m = 1, K = 2, L = 2,3, ... gamma band, 
n = 1, n' = 0, m = 0,K = 0,L = 0,2,... beta band. (3.12) 

In the situation when the /3 potential is an infinite square well, the wave function has the 
expression: 



\ns;n'm) = fUP)Sn'Ucos3^ y {D^jj, + (-l)^+^/^^,,_^) ,m = -. (3.13) 

with f^r given by Eq.(2.10) and the irrational index s given by Eq. (13.61) and S'„/m(cos37) 
defined by Eq.(2.39). The quantum numbers defining the ground, beta and gamma bands 
are as follows: 

^ = 1, n = 0, m = 0, K = 0, L = 0,2, ... ground band, 

^ = 1, ,n' = l, m = 1, K = 2, L = 2,3, ... gamma band, 

^ = 2, ,n' = 0, m = 0,K = 0,L = 0,2,... beta band. (3.14) 
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Once the wave functions are determined by solving the corresponding eigenvalue equa- 
tions, we can proceed to calculating the electric transition probabilities. In order to get a 
feeling about how sensitive the matrix elements of 7 depending terms of the transition op- 
erator are to changing the parameter c, we have plotted them in Figs. 10, 11, versus c. From 
there one notices that in a large interval of c, the matrix elements are slowly varying with 
c. The diagonal matrix elements of cos 7 (the first panel in the left column) are shghtly in- 
creasing by 0.01 starting with the values 0.823, 0.844 and 0.854 at c = 0. The corresponding 
matrix elements of sin 7 are changing just a little when we vary c, starting with the values 
0.475 0.487 and 0.493. The magnitudes of the matrix elements between states belonging to 
the same multiplet (see Fig.l) are small. The matrix elements characterized by the same 
An = 1 are relatively large for cos 7 but small for sin 7. 

The reduced E2 transition probabilities have been calculated by using alternatively a 
harmonic, T^'*\ and an anharmonic transition operator, T2^^\ having the expressions: 

7^;) = tfi [cos^Dl, + '^iD% + Dl_,)j , 

hf^P' (- ^«^27i^2o + ^(^i + ^^,-2)) (3.15) 

The strengths t, ti and ^2 are free parameters which are fixed by fitting one and two particular 
B{E2) values, respectively. Due to the structure of the wave functions specified above, 
the matrix elements between the states involved in a given transition are factorized into 
matrix elements of the transition operators factors depending on /3, 7 and the Euler angles, 
respectively. 
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FIG. 6: For three sets of parameters {ui,U2) = 
(-1, -1); (100, 1100); (5, 85)keV , the 7 potential 
[7(7) (3.8), with (/3^) = 1, is plotted as a func- 
tion of 7 in the pannels a), b), c), respectively. 



FIG. 7: The 7 potentials (3.8), cor- 
responding to the parameters (^1,^2) = 
(10, -20); (-5, -85); (-20, -40)keV and {P^) = 
1, are represented as functions of 7 in panels d), 
e) and f), respectively. 
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FIG. 8: The spheroidal functions corresponding 
to the 7 potentials from Figs. 6 a), 6 b) and 6 
c), are given as functions of 7 in panels a), b) 
and c) respectively. The functions describe the 
ground band states of angular momenta 0, 4 and 
8, respectively. 



FIG. 9: The spheroidal functions corresponding 
to the 7 potentials from Figs. 6 d), 6 e) and 6 
f), are plotted versus 7 in panels d), e) and f) 
respectively, for three states of angular momenta 
0, 4 and 8 belonging to the ground band. 
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FIG. 10: The matrix elements of cos 7 
{Tnl.n'm') ^nd sin 7 are plot- 

ted as functions of the parameter c in- 
volved in Eq.(2.39), for {nm;n'm') = 

(11; 11), (00; 00), (10; 10) (1st and 3rd pan- 
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els) (GO; 10), (00, 11), (11, 10) (2nd and 4th 



FIG. 11: The same as in Fig. (10) 
but for different quantum numbers 
{nm;n'm') = (30; 20), (31; 21) (1st and 3rd 
panels) (33, 32), (21; 20), (22; 21), (32; 31) (2nd 
and 4th panels). 



IV. NUMERICAL RESULTS 



The formalism described in the previous Section has been apphed for ^^°Nd, ^^^Gd and 
^^^Os. The choice is justified by the values of the ratio of the excitation energies for the first 
two excited states in the ground band. Indeed, these are 2.93, 3.015 and 2.82 respectively 
and therefore they are expected to have the features of X(5) symmetry. As we stated in 
Introduction the present formalism is close to the X(5) symmetry. Indeed, it goes to X{5) in 
the limit 7^0 and, on the other hand, the spheroidal function equation has been derived 
by expanding the gamma depending terms in the initial Hamiltonian in terms of sin 37. 
However, it brings two new things namely i) the wave function is periodic in 7, the matrix 
elements of 7 depending functions being performed with the measure | sin 37 1^7 and ii) the 
factor function describing the degree of freedom satisfies a Schrodingcr equation involving 
the Davidson's potential. Numerical calculations will show us what arc the quantitative 
corrections to the X{5) picture brought by curing the drawbacks of the preceding approaches. 

The 7 potentials for the three nuclei have been chosen from those given in Fig. 6. Indeed, 
the potentials ^^^Gd, ^^°Nd and ^^^Os are those from Figs. 6 a), c) and b) respectively. The 
spheroidal functions corresponding to these potentials are represented in Figs 8 a), 8 c) and 
8 b), respectively. We varied the shapes of gamma potential and made the option for that 
one which yields a good agreement between the calculated B(E2) values associated to the 
transitions from gamma to ground band and the corresponding experimental data. 

The energy levels in the three bands are given in units of E2+ and therefore we need only 
the ratio of the parameters A and F involved in the energy expression given by Eq. (3.10). 
The term Eq is not depending on the quantum number defining the states but it is considered 
to depend on band. In our calculation Eq for ground and beta bands are taken to be equal 
and therefore they do not affect the relative energies in the two bands, while for gamma 
band was fixed so that the head state energy is recovered. The parameters defining the 
transition operator have been fixed by fitting the B(E2) values for the transitions 2+ 0+ 
and 2+ — > 0^. The results obtained in this manner are collected in Table I. 

Now let us proceed at describing separately the results for each of the three nuclei con- 
sidered here. 

The calculated energies for ^^'^Nd are listed in Table II. The results of present paper are 
given in the columns D and ISW. They have been obtained by using a spheroidal function 



30 





isoNd 


154q^ 


1920s 


ui[keV] 

-L L J 


100 


-1 


5 


U2[keV] 


1100 


-1 


85 


Eo\keV] 


-27.696 


-222.243 


9.635 




-3.545 


-12.892 


0.367 


F/A 


0.532 


17.728 


0.191-10-3 




0.0328 


1.027 


2.896-10-3 




147.889 


124.065 


90.583 




161.909 


135.831 


99.174 




209.443 


184.933 


120.976 




228.706 


219.473 


138.962 




100.178 


106.240 


50.077 




116.624 


156.629 


71.225 



TABLE I: The parameters ui,U2,Eo,F/A, involved in the energy expression (3.10), calculated 
by the method described in the text, are given for i^'^Nd, ^^^Gd, ^^^Qg^ Also we give the values 
for the parameters t and ti,t2 involved in the harmonic and anharmonic transition operators, 
respectively. They were obtained by fitting the B(E2) values for the transition 2^ 0^ if the 
harmonic transition operator is used, and for the transitions 2^ 0^ and 2+ 0^ for the 
anharmonic transition operator. 



description in the 7 variable while for the f3 potential, an infinite square well (ISW) and the 
Davidson's potentials (D), have been alternatively used. These results are compared with 
the corresponding experimental data {47! . 48], listed in the column headed by Exp., as well 
as with the theoretical results obtained within the X{5) approach. If we enlarge the list 
for each band the deviations of predictions from the corresponding experimental data are 
increasing functions of angular momentum. 
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State 


Exp. X(5) ISW D 


2+ 


1111 




2.93 2.90 2.93 2.93 




5.53 5.43 5.53 5.53 




8.68 8.48 8.70 8.73 


10+ 


12.28 12.03 12.42 12.50 




16.27 16.04 16.66 16.82 




5.19 5.65 5.30 4.20 




6.53 7.45 7.05 5.45 




8.74 10.69 10.24 7.65 


6+ 


11.84 14.75 14.27 10.64 




8.16 8.16 8.16 8.16 


3+ 


9.22 8.93 9.03 9.35 




10.39 9.83 10.08 10.70 



TABLE II: Excitation energies of some states from ground, beta and gamma bands of ^^'^Nd, 
given in units of the 2+ energy, obtained with three different approaches, -^(5), infinite square 
weh potential (ISW) for /3 and spheroidal functions formalism for 7, Davidson's (3 potential plus 
spheroidal functions method for 7 (D), are compared with the corresponding experimental data. 



The intraband ground band and beta band transitions as well the gamma to ground and 
beta to ground transitions were calculated by using both a harmonic and an anharmonic 
structure for the transition operators with the strength parameters from Table I. The final 
results are those from Table III. In the mentioned Table we give also the corresponding 



experimental data, taken from Ref. 



47 



48j, as well as the results provided by the X(5) 
formalism. It is interesting to remark that except for the gamma to ground transitions, 
the harmonic approach of the ISW calculations provides identical results with the X(5) 
calculations. We may say that the agreement with the data is improved by adding the 
anharmonic effects for some transitions but for some other the discrepancies are increased. 
In an overall analysis, the agreement is improved by anharmonicities. 
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B{E2-J+ 




Exp. X(5) 


ISW 


D 






115 


115 


115 


115 


115 


115 


4+ ^ 




182 


184 


184 


177 


197 


182 






210 


228 


228 


210 


266 


226 






278 


262 


262 


233 


334 


260 


10+- 




204 


288 


289 


249 


394 


387 






114 


92 


92 


91 


148 


121 




2+ 


170 


138 


138 


136 


210 


167 






39 


72 


72 


43 


105 


58 






1.2 


2.4 


2.4 


0.4 


0.6 


0.1 


2+ ^ 


2+ 


9 


10 


10 


4 


9 


3 


2+^ 


4+ 


17 


42 


42 


26 


69 


36 


4- 


2+ 


0.12 


0.56 


1 


0.01 0.61 


4.6 






7 


7 


7 


3 


4 


0.3 


4- 


n 


70 


32 


32 


18 


51 


22 


2+ ^ 




3.0 


3.4 


0.6 


3 


0.6 


3 


2+ ^ 


2+ 


5.4 


5.1 


0.9 


4.7 


1 


5.0 


2+ ^ 

^7 


4+ 

9 


2.6 


0.3 


0.05 


0.3 


0.05 


0.3 


4+ ^ 

7 


2+ 


0.9 


2.3 


0.4 


2.1 


0.4 


2.3 


4+ ^ 

7 


4+ 

9 


3.9 


7.3 


1.3 


6.8 


1.5 


8.2 



TABLE III: The B(E2) values for i^^Nd calculated in three different formalisms, X(5), ISW, B, and 
the corresponding experimental data are given in units of e^/m^ x 10^. The results of the present 
work are obtained by using spheroidal functions for 7 and alternatively an infinite square well 
potential {ISW) and Davidson's potential (D) for the variable /3. The results for X(5) formalism 
are taken from Ref . . In the first columns headed by ISW and D respectively, are the results 
obtained with a harmonic quadrupole transition operator. In the second columns ISW and D we 
give the results obtained with an anharmonic quadrupole transition operator. Calculations made 
with the X{5) formalism correspond to a harmonic transition operator. 
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State 


Exp. X(5) ISW 


D 


CSM CSM2 


2+ 


1 1 1 


1 


1 1 




3.015 2.90 3.015 3.015 


3.110 2.932 




5.83 5.43 5.84 


5.83 


6.105 5.612 




9.30 8.48 9.39 


9.35 


9.835 8.936 


10+ 


13.30 12.03 13.55 13.81 


14.188 12.850 




5.53 5.65 4.12 


3.35 


5.662 5.335 




6.63 7.45 5.71 


4.51 


6.413 6.025 




8.51 10.69 8.70 


6.68 


8.101 7.714 




11.10 14.75 12.63 


9.75 


10.626 10.325 




14.27 19.44 17.36 13.58 13.879 13.756 




17.83 24.69 22.79 18.15 17.782 17.916 


2+ 


8.10 8.10 8.10 


8.10 


7.771 7.817 




9.16 8.86 9.00 


8.89 


9.047 8.528 




10.27 9.75 9.87 


10.33 10.073 9.514 


5+ 


11.64 10.75 11.01 


11.74 11.301 10.604 


6+ 


13.05 11.86 12.31 


13.31 


12.778 12.002 


7| 


14.71 13.07 13.75 15.02 14.386 13.421 



TABLE IV: The same as in Table II but for ^^^Gd. 
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Although ^^^Gd is a deformed nucleus, having the quadrupole deformation ^ 0.25 [56|, 
the authors of Ref . {49 1 consider it as a good candidate for the critical point in a phase 
transition which takes place along the chain of the Gd isotopes. This shape transition has 
been studied recently by two of us (A.A.R. and A. F.) in Ref. 53| within the Coherent State 
Model (GSM). Here the results obtained through ISW and D formalism are compared to 
those obtained within other phenomenological approaches like X(5), GSM, and GSM2. The 
GSM2 differs from GSM by the model states for the beta band. A full list of references 
concerning GSM may be found in Ref.j^. The results for X(5) are taken from Ref. 57] with 
a suitable scaling when one passes from ^^^Sm to ^^^Gd. We notice that the D formalism 



provides energies which are closest to the experimental data 



51 



52|. Goncerning the 



B(E2) values, from Table V one notices that D and ISW formalisms are reproducing better 

nn 

the data [50, 151| for the intraband ground and beta bands transitions. As regards the beta 
to ground transitions the X(5) formalism yields an overall better agreement with the data. 
Also, we note that the approach D is describing quite well the transitions (J + 2)^ but 

its predictions for are smaller than the corresponding data by a factor of about 

10. By contrary the theoretical transitions ^ ( J — 2)^ are larger by a factor 10 to 50 
than the experimental results. The best description for the interband transitions collected in 
Table V is provided by the ISW approach. The anharmonic term of the transition operator 
brings important contribution to both the intra and interband transitions. 
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Transition Exp. X(5) 


ISW 


D CSM 


CSM2 


2+ 






77.3 


77.3 


77.3 


77.3 


77.3 


77.3 77.3 


77.3 


77.3 77.3 


4+ 

9 






117.8 123.5 123.7 117.8 132.4 117.8 114.4 114.9 117.2 117.3 








138.2 


153 


153.4 138.9 179.1 


138.9 133.1 


133.9 139.7 139.9 








152.6 176.1 


176.3 152.8 


225 


152.2 148.7 148.4 159.4 159.8 


10+ 




► 8+ 173.1 193.8 194.6 162.4 265.6 


160 163.7 163.8 178.4 178.9 








49 


61.7 


61.5 


61.3 


99.2 


72.7 103.8 


67.1 


102.8 103.8 


4 






122 


92.9 


92.8 


91.2 


141.1 


99.3 152.4 


99.5 


150.9 152.4 






4 


111 


113.6 113.6 109.5 152.4 107.5 175.6 116.2 


174 175.6 






2+ 


25.8 


48.8 


48.2 


26 


70.6 


26.5 l(-3) 


1.32 


19.9 20.0 


2+ 




2+ 


4.0 


7.0 


6.4 


2.1 


6.1 


0.44 2(-4) 


0.34 


4.13 4.16 






4+ 


11.9 


27.9 


28.3 


15.4 


46.4 


15.9 3(-3) 


0.76 


12.98 13.08 








0.35 


0.54 


0.73 


0.06 


0.4 


5.4 


0.21 


2.17 2.19 








3.8 


4.83 


4.72 


1.29 


2.8 


0.044 3(-3) 


0.33 


3.44 3.47 








12 


21.5 


21.6 


10.7 


34.3 


8.16 8(-3) 


0.68 


13.42 13.53 


6/^ 






0.27 


0.50 


0.50 


0.18 


7 


17.3 l(-4) 


0.17 


1.70 1.71 



TABLE V: The same as in Table III but for ^^^Gd. The notation k{—m) stands for the number 
k ■ 10-"^. 
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Li, L 




Exp. 


X(5) 


ISW 


D 


CSM 


CSM2 


% 


2^ 


2^ 




1.52 


1.60 


1.60 


1.52 


1.71 


1.52 


1.48 


1.48 


1.52 1.52 






49 


29 


1.17 


1.24 


1.24 


1.18 


1.35 


1.18 


1.16 


1.16 


1.19 1.19 


89 




6. 


4. 


1.10 


1.15 


1.15 


1.10 


1.26 


1.10 


1.12 


1.11 


1.14 1.14 


10, 




8. 


6. 


1.13 


1.10 


1.10 


1.06 


1.18 


1.06 


1.10 


1.10 


1.119 1.113 


27 


% 


27 


2. 


0.468 


0.666 


0.654 


0.642 


0.618 


0.594 


0.509 


0.468 


0.501 0.468 


2, 


4. 


2, 


29 


0.144 


0.052 


0.053 


0.054 


0.054 


0.057 


0.087 


0.131 


0.068 0.069 


3'-y 


2. 




4^ 


1.006 


2.368 


2.335 


2.274 


2.208 


2.096 


1.302 


0.975 


1.432 1.289 




2. 


47 


4^ 


0.148 


0.315 


0.311 


0.304 


0.286 


0.273 


0.159 


0.126 


0.152 0.123 


47 


6. 


47 


4. 


0.27 


0.088 


0.090 


0.092 


0.092 


0.097 


0.377 


0.377 


0.115 0.117 




4. 




6. 


0.744 


1.667 


1.655 


1.619 


1.560 


1.486 


0.657 


0.401 


0.786 0.665 


6^ 


4. 


6^ 


% 


0.081 


2.333 


0.25 


0.246 


0.229 


0.219 


0.081 


0.046 


0.076 0.050 


2„ 


2/3 


27 


2, 


1.00 


0.032 


0.054 


0.108 


0.061 


0.145 


1.206 


0.751 


2.238 0.048 


2/3, 


0. 


2^ 


2^ 


0.123 


0.429 


0.257 


0.067 


0.069 


1.716 


0.0 


0.561 


0.538 0.475 


2^, 


4. 


2/3 


2. 


2.76 


4.00 


4.45 


7.20 


7.6 


36 


13.24 


2.257 


3.141 3.141 


4/3, 


2. 


4/3 


4. 


0.086 


0.112 


0.155 


0.046 


0.145 


121 


0.001 


0.625 


0.630 0.630 


4/3 


6. 


4/3 


4^ 


2.63 


4.45 


4.57 


8.30 


12 


184 


2.484 


2.073 


3.896 3.896 


6/3, 


4. 


6/3 


6p 


0.08 




0.13 


0.22 


2475 


292 


0.0 


0.071 


0.555 0.555 


2/3, 


Op; 


2^, 


0/3 


0.008 


0.049 


0.027 


0.002 


0.004 


0.568 


0.0 


0.0028 0.022 0.019 


4/3, 


2.; 


4^, 


2/3 


0.0025 0.0058 0.0079 0.0006 0.0029 0.530 


0.0 


0.0021 0.014 0.014 




4.; 


6/3, 


4/3 


0.0024 0.0044 0.0044 0.0017 


0.045 


0.676 


0.0 


0.0015 0.010 0.010 


8/3, 


6^; 


8/3, 


6/3 


0.006 




0.003 


0.002 


0.204 


2.147 0.977 0.0011 0.007 0.007 



TABLE VI: Calculated branching ratios B{E2;Lf Lp/B{E2; L^+) (denoted by 

Li,Lf;L[,Lf' ) for some interband as well intraband transitions in ^^^Gd are compared with the 
corresponding experimental data. The results of this paper, labeled by ISW and D, were obtained 
by using spheroidal functions for 7 and for the /3 variable an infinite square well and the Davidson's 
potentials, respectively. Results from the first columns headed by ISW and D were obtained by 
using a harmonic transition quadrupole operator while in the second columns are listed the results 
corresponding to an anharmonic transition operator. 
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Starea Exp. 


X(5) 


ISW 


D 


2+ 


1 


1 


1 


1 


4+ 

9 


2.82 


2.90 


2.90 


2.90 




5.29 


5.43 


5.43 


5.43 




8.30 


8.48 


8.48 


8.48 


10+ 


11.75 12.03 12.03 12.01 




15.60 16.04 16.04 16.01 




4.65 


5.65 


5.65 


4.44 


2+ 


5.48 


7.45 


7.45 


5.73 


2+ 


2.38 


2.38 


2.38 


2.38 




3.35 


2.60 


3.24 


3.32 




4.42 


2.87 


4.28 


4.44 


5+ 


5.56 


3.16 


5.48 


5.75 


6+ 

7 


7.12 


3.48 


6.81 


7.20 




8.32 


3.84 


8.27 


8.81 


8+ 


10.37 


4.22 


9.86 


10.54 


10+ 


14.06 


5.07 


13.41 


14.40 



TABLE VII: The same as in Table II but for ^^^Os 

Some branching ratios of the intraground band as well as of the interband transitions 
are given in Table VI. The failure of the D approach to describe these branchings for the 
P ^ g transitions is noticeable. The reason consists in the fact that within this formalism the 
predictions for the transitions Jg — > J+ are too small comparing them with the experimental 
data. The other ratios are described reasonable well by all theoretical models. 

The results for ^^^Os are presented in Tables VII and VIII. Therein, we give also the 
experimental data and the theoretical results yielded by the X{5) formalism. First we note 
two specific features for this nucleus: i) E4+/E2+ — 2.82, which recommends it as a good 
candidate for the ^(5) symmetry and ii) £'2+ < Eq+. This inequahty characterizes the 
gamma unstable nuclei. Prom Table VII we remark the very good description of energy 
levels by the D formahsm. Also, from there it results that the ^(5) formahsm fails to 
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Transition 


Exp. 


X(5) ISW D 


2+ 






0.424 


0.424 0.424 0.424 0.424 0.424 


4+ 

9 






0.497 


0.678 0.678 0.656 0.726 0.673 








0.660 


0.841 0.840 0.787 0.981 0.836 








0.754 


0.966 0.964 0.879 1.23 0.966 


10+ 






0.688 


1.060 1.06 0.947 1.45 1.06 


4+ 




2+ 

-^7 


0.298 


0.269 0.275 0.274 0.302 0.288 


6+ 




47^ 


0.336 


0.611 0.610 0.590 0.730 0.643 


8+ 




6+ 


0.314 


0.823 0.807 0.761 1.05 0.858 


2+ 

^7 






0.037 


0.012 0.009 0.037 0.009 0.037 


2+ 

^7 




2+ 


0.303 


0.019 0.014 0.057 0.015 0.062 


4+ 




2+ 

7 


0.014 


0.001 0.006 0.026 0.007 0.029 


4+ 

7 






0.002 


0.008 0.006 0.026 0.007 0.028 


^7^ 






0.203 


0.027 0.020 0.084 0.023 0.102 








0.012 


0.006 0.006 0.025 0.007 0.032 


6+ 

"7 






0.0004 0.006 0.006 0.025 0.007 0.031 


6+ 

"7 




6.^ 


0.171 


0.023 0.023 0.100 0.029 0.138 



TABLE VIII: The same as in Table III but for ^'^^Os. The units for the B(E2) values are e^fe^. 

describe the excitation energies in the gamma band. Concerning the E2 transitions the data 
from Table VIII show that the ISW formalism with a harmonic transition operator yields 
results very close to the ones produced with the ^(5) formalism. However including an 
anharmonic term in the expression of the transition operator one obtains a slightly better 
agreement with the experimental data. 

Concluding the numerical analysis of this Section wc may say that the formalism presented 
here has not only the merit of removing two drawbacks of the previous X(5) model (the 
functions in 7 are not periodic while the Hamiltonian is not Hermitian.) but also provides 
a better quantitative description. 
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V. CONCLUSIONS 



Here we shall summarize the main results presented in the present paper. Starting with 
the differential equation for (5 and 7 variables, involving also the Euler angles, provided by 
the liquid drop model, one may define analytically solvable equations for both deformation 
variables. This is possible under certain circumstances which are described in details in 
Sections II and III. The first model which achieved this situation corresponds to the so 
called -^(5) symmetry and has the merit of describing in a very simple fashion the properties 
of the critical point in shape phase transition. Two specific features are considered as 
drawbacks of the model: l)the function describing the variable 7 are non-periodic functions 
and moreover are normalized to unity on a non-bounded interval. These two properties 
confiict the symmetry properties required by the starting liquid drop Hamiltonian. 2)In 
particular if we preserve the integration measure for 7 to be | sin 2>'^\d'^ then the Hamiltonian 
depending on 7 is not Hermitian. 

The scope of this paper was to remove these two drawbacks. The solution offered here is 
to use in the 7 space a Hamiltonian which has the spheroidal functions as eigenstates. These 
functions are, indeed, periodic in the interval [0,27r]. Moreover, the model Hamiltonian is 
Hermitian with respect to the integration measure | sin 87! 0(7. Another solution would be the 
use of a 7- Hamiltonian which admits the Mathiew functions as eigenfunctions. However, 
applications based on this solution are postponed for another publication. We have proved 
that both solutions lead the X(5) model in the limit of I7I small. As regards the (3 variable, 
the description is performed with the generalized Laguerre functions which are solutions 
of a Schrodinger equation involving the Davidson's potential. A particular case of this 
description is the situation when the centrifugal terms is vanishing, i.e. (5q — 0. Of course, 
that is the oscillator potential in (3. 

It is inferred that each of the solvable models in the j3 variable, which have been previously 
used as E{h) models, may be associated to the spheroidal function description by separating 
the term proportional in and not depending on 7 to renormalize the equation in (5. Two 
examples are given in this paper where the Davidson potential and an infinite square well 
potential are used. Associations of other (3 potentials and the spheroidal function approach 
arc presently under our consideration. 

The numerical applications to ^^°Nd, ^^''Gd and ^^^Os reveal a good agreement with 
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experimental data. Moreover, the comparison with the results yielded by X(5) calculations 
suggest that the present approach provides a slightly better quantitative description of the 
data. For most of the data presented here the results of the present approach are close 
to those obtained with the X(5) formalism. However, this is not a surprising feature if 
we recall that the 7-Hamiltonian has been derived through a sin3'j expansion approach. 
If we consider higher angular momentum states in the three major bands, the deviations 
of energies predicted by the X(5) calculations from the corresponding experimental data 
are larger than those obtained with the present formalism. One notes that the gamma 
band energies are better described by the present approach. As a matter of fact the X(5) 
calculations fail to reproduce the excitation energies in the gamma band of ^^^Os. In the 
case of ^^"^Gd wc compared the theoretical results of this paper with those obtained with two 
versions of the Coherent State Model. The qualities of the agreement with the experimental 
data obtained with the three sets of calculations are comparable with each other. The 
essential difference is that CSM works quite well for all even Gd isotopes while the -'^(S)- 
type models work only for the critical point of the shape phase transition. 

Perhaps this is the prize we have to pay by endorsing the variable separability. The model 
beauty cannot substitute the virtue of the variable mixing to account for details specific for 
one phase or another. 
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VI. APPENDIX A 



In order to help the reader to check the expressions given in the text, here we give some 
intermediate results concerning the partial expansions in 7, used in deriving the fourth order 
expansion for the Hamiltonian considered in Section II. Thus, the useful expansions are listed 
below: 

1 4 87 872 I673 1047^ ,5 



sin2(7-f) 3 3^3 3 3v^ 27 

1 _ 4 , 87 , 87' , 16o'-^ , 10474 



sin2(7+f) 3 3^3 3 3^3 27 
1 1 1 7^ 27'^ 



S = ^ + 3 + 15 + 189+^^'' 



sm 7 7 

^ = A + i + ^ + ^ + ObP. (3.1) 

sm (37) 97^ 3 5 7 

In order to perform the expansion around 7r/6 one needs the following expansions in terms 
of y = I7 — 7r/6|: 



sin2(7 - f ) ' " ' 3 



1 



r- o 176t/3 728|/4 

4 + 8^32/ + 402/2 + + — ^ + 0[2/]^ 



sin2(7 + f) ' V ^ , ^ , ^ , g 



sm^7 V3 3 

1 l + 9y2 + 54y4 + e»[y]5. 



sin (37) 



(3.2) 



VII. APPENDIX B 

In the rotational term: 



we shall consider the identity: 



1^1 9 

4 S sin' (7 - f A;) " sin2(37)- ^^"^^ 
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The term W{''y, Q) acquires a more convenient form: 

9 1 



sin^(37) sin^7 
1 



iQi + Qi + Qi) + 

1 



+ 



sin (87) sin 7 



Ql 



_sin^(7-f) sin^(7-f,j 
In the regime of I7I -C 1 one uses the expansions 

9 1 3 2772 

= — = \ u 

472 ^ 4 



{Qi - Qi 



(B.3) 



4 sin 37 
1 



+ . + ^ + 0(7') 



sin^ 7 



cos 37 



1 1 

7 

1 - ^ + Oi^') 



COS237 = 1-972 + 0(7^), 
in connection with the expression IB.3I of W{'-f,Q). The result is: 

+ ^{Ql-Qi) + o{i'). 



(B.4) 



Inserting this expression into the Hamiltonian 



H-E 



^ 9 
972 ^ 4 



1 + 



sin 37 



?7(7)-lV(7,g), 



(B.5) 



(B.6) 



and averaging the result with an axially symmetric rotational state, for which (Qi) = (g2), 
and denoting the result by Hq, one obtains: 



92 9 

+ T-^(7)-n7), 



F(7) 



8sin2(37) 



\ sin^ 7 



L(L + 1) - 3g^ 



(B.7) 
(B.8) 



972 4 
'l{L + 1) - - 2 

Here we used the notation = (Qk)- If the averaging is performed with a Wigner function 
D^jj^, the result for V{'y) would be: 

9 1 



L{L + l)-2-K' 



(B.9) 



4sin^7^ ' ' 8 [sin2(37) sin2 7 

In the limit of I7I <^ 1 the above expression of V{'y) can be expanded in powers of 7. The 
second order expansion is: 



12 V37^ 



+ 1 + W]{K^-1) + - L{L + 1] 



(1 + 27^) + 0(73). (B.IO) 
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On the other hand making use of the expansion 

9 1 



.20 - . 2 +;T(l + 2sin2 7) + 0(73), 
sm 37 sm 7 3 



one obtains the following second order expansion in sin 7: 

4 sm 7 3 

From this expression one obtains immediately the expansion in sin (37) 



L{L + l)-K^-2 (1 + 2 sin^ 7) + Oi-f^). 



+ 1) -K^ -2 



l + -sin2 37) +0(7' 



(B.ll) 



(B.12) 



[B.13) 



4 sin^ 37 3 

The three expansions for V{'y) are useful to study different representations for the wave 
function in 7. Thus, inserting (IB.lOp in Eq.(B.7) one obtains a differential equation for 
the Laguerre functions. Here, all corrections in 7^ were included. Using the expansion in 
sin 7 the differential equation for the variable 7 becomes an equation for the spheroidal 
function after the change of variable x = cos'y is performed. Finally, we mention that the 
use of Eq. flB.13P leads to a differential equation for the spheroidal functions in the variable 
X = cos(37). 



VIII. APPENDIX C 



The renormalization of the equation in j3 due to the terms coming form the rotational 
term, is based on the expansion in sin 37 given by Eq. (B.13). Multiplying 1^(7) given by 
the quoted equation with one obtains terms which depend on 7 and terms which do not 
depend on this variable. The latter terms are: 

1 



R 



3/32 



Our option for renormalization is based on the approximation: 



R 



\LiL + l)- ^ 



(C.l) 



(C.2) 



3/52 ^ ' 3(/52) 

In this way the indices for the generalized Laguerre functions, if we use the Davidson po- 
tential, or of Bessel function if we use an infinite square well 7 potential are those used 
in the present paper and Ref. 59|]. If the separation of the two types of terms is achieved 
differently: 



R 



3/3^ 



2] 



3(/32) 



(C.3) 
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then the irrational indices for Laguerre and Bessel functions are: 



3 1 / 3 



Note that in the case the /3 potential is a harmonic oscillator and the 7 potential is ignored, 
the separation of the (3 variable takes place in a natural manner, the 7 and Euler angles 
depending terms being just the Casimir operator of the group R{5). Due to this separation 
the f3 wave function is not depending on the quantum number K. Here such a simple 
picture does not hold any longer and a K depending term may renormalize the equation for 
13. Therefore, it becomes meaningful to consider the full term R (IC.ll) for the renormalization 
purpose. In this case the two irrational indices for the generalized Laguerre function and 
the Bessel function are: 



2 V 3 V 2/ 3 2 
2 + V3(^+2) -^+2- (^-^^ 



The difference between the index s given by Eq. (IC.5I) and that used in Ref . 57| is caused by 
the fact the therein the term from Eq. OC.ip is not used for renormalizing the equation 
for f3. 



[1] A. Bohr, Mat. Fys. Medd. Dan. Vid. Selsk. 26 (1952) no.l4; A.Bohr and B.Mottelson, Mat. 

Fys. Medd. Dan. Vid. Selsk. 27 (1953) no. 16 
[2] A. Faessler and W. Greiner, Z. Phys. 168 (1962) 425; 170 (1962) 105; 177 (1964) 190; A. 

Faessler, W. Greiner and R. Sheline, Nucl. Phys. 70 (1965) 33. 
[3] G. Gneus, U. Mosel and W. Greiner, Phys. Lett. 30 B (1969) 397. 

[4] P. Hess, J. Maruhn and W. Greiner, Phys. Rev. C23 (1981) 2335; J. Phys. G, 7 (1981) 737. 
[5] A. A. Raduta, V. Ceausescu, A. Gheorghe and R. M. Dreizler, Phys. Lett. 99B (1981) 444; 

Nucll. Phys. A381 (1982) 253. 
[6] A. A. Raduta, V. Ceausescu and A. Faessler, Phys. Rev. C36 (1987) 2111. 
[7] A. A. Raduta, C. Lima and A. Faessler, Z. Phys. A 313 (1983) 69. 

45 



[8] A. A. Raduta, Al. H. Raduta and A. Faessler, Phys. Rev. C55 (1997) 1747; A. A. Raduta, D. 

lonescu and A. Faessler, Phys. Rev. C65 (2002) 064322. 
[9] A. A. Raduta, in Recent Res. Devel. Nuclear Phys.,1 (2004):l-70, ISBN:81-7895-124-X. 
[10] L. Wilets and M. Jean, Phys. Rev. 102 (1956) 788. 
[11] A. S. Davydov and G. F. Filippov, Nucl. Phys. 8 (1958) 788. 
[12] A. Arima and F. lachello, Ann. Phys.(N.Y.) 99 (1976) 253; 123 (1979) 468. 
[13] F. lachello and A. Arima, The Interacting Boson Model (Cambridge University Press, Cam- 
bridge, England, 1987). 

[14] R. F. Casten, in Interacting Bose-Fermi Systems in Nuclei, edited by F. lachello (Plenum, 

New York, 1981), p. 1. 
[17] D. J. Rowe, Nucl. Phys. A 745 (2004) 47. 
[16] P. S. Turner, D. J. Rowe, Nucl. Phys. A 756 (2005) 333. 
[17] G. Rosensteel, D. J. Rowe, Nucl. Phys. A 759 (2005) 92. 
[18] J. H. Ginocchio and M. W. Kirson, Phys. Rev. Lett.44(1980) 1744. 
[19] A. E. L. Dieperink, O. Scholten and F. lachello, Phys. Rev. Lett.44 (1980) 1747. 
[20] F. lachello, Phys. Rev. Lett. 85 (2000) 3580. 
[21] F. lachello, Phys. Rev. Lett. 87 (2001) 052502. 
[22] R. F. Casten and N. V. Zamfir, Phys. Rev. Lett. 85 (2000) 3584. 
[23] R. F. Casten and N. V. Zamfir, Phys. Rev. Lett. 87 (2001) 052503. 

[24] R. M. Clark, M. Cromaz, M. A. Deleplanque, M. Descovich, R. M. Diamond, P. Fallon, I. 

Y. Lee, A. O. Macchiavelli, H. Mahmud, E. Rodriguez- Vieitez, F. S. Stephens and D. Ward, 

Phys. Rev. C 69 (2004) 064322. 
[25] N. V. Zamfir et al., Phys. Rev. C65 (2002) 044325. 
[26] Da-U Zhang and Yu-xin Liu, Phys. Rev. C65 (2002) 057301. 
[27] L. Fortunato and A. Vitturi, Jour. Phys. G:Nucl. Part. Phys. 29 (2003) 1341. 
[28] Dennis Bonatsos, D. Lenis, N. Minkov, D. Petrellis, P. P. Raychev, P. A. Terziev, Phys. Lett. 

B 584 (2004) 40. 
[29] P. M. Davidson, Proc. R. Soc. 135 (1932) 459. 

[30] A. A. Raduta, A. Gheorghe and A. Faessler, J. Phys. G:Nucl. Part. Phys.,31 (2005) 337. 

[31] A. A. Raduta, F. D. Aaron and I. I. Ursu, Nucl. Phys. A 772 (2006) 20. 

[32] M. A. Caprio, P. Cejnar and F. lachello. [arXiv:0707.0325V l [quant-ph] . Annals of Physics NY, 

46 



323 (2008) 1106. 

[33] L. Fortunate, |arXiv:nucl-th/0411087|vl; Eur. Phys. J. A 26 (2005) sOl. 

[34] A. Gheorghe, A. A. Raduta and V. Ceausescu, Nucl. Phys. A 296 (1978) 228. 

[35] A. A. Raduta, A. Gheorghe and V. Ceausescu, Nucl. Phys. A311 (1978) 118. 

[36] A. Gheorghe, A. A. Raduta and V. Ceausescu, Nucl. Phys. A637 (1998) 201. 

[37] E. Chacon, M. Moshynsky and R. T. Sharp, J. Math. Phys. 17 (1976) 668. 

[38] T. M. Corrigan, F. J. Margetan and S. A. Wilhams, Phys. Rev. C 4 (1976) 2279. 

[44] D. J. Rowe, Nucl. Phys. A 735 (2004) 372. 

[40] D. Bes, Nucl. Phys. 210 (1959) 373. 

[41] A. Gheorghe, A. A. Raduta and A. Faessler, Phys. Lett. B 648, 171, (2007). 

[42] P. M. Davidson, Proc. R. Soc. 135 (1932) 459. 

[43] J. P. Elliot, J. A. Evans, P. Park, Phys. Lett. B 169 (1986) 309. 

[44] D. J. Rowe, C. Bahri, J. Phys. A 31 (1998) 4947 

[45] M. Abramowitz, and I. A. Stegun, (Eds.), Handbook of Mathematical Functions with Formu- 
las, Graphs, and Mathematical Tables, 9th printing. New York: Dover, 1972, pp. 751-759. 

[46] S. De Baerdenmacker, L. Fortunato, V. Hellemans, K. Heyde, Nucl. Phys. A 769 (2006) 16. 

[47] R. Kriicken et al. Phys. Rev. Lett. 88 (2002) 232501. 

[48] E. der Mateosian, J. K. Tuh, NDS 75 (1995) 827. 

[49] R. M. Clark et al, Phys. Rev. C68 (2003) 037301. 

[50] P. O. Lipas et al Phys. Scr. 7 (1983) 8. 

[51] C. Girit, W. D. Hamilton and C. A. Katfas, J. Phys. G: Nucl. Part. Phys.9 (1983) 797. 

[52] C. W. Reich, R. G. Helmer, NDS 85 (1998) 171. 

[53] A. A. Raduta and Amand Faessler, J. Phys. G: Nucl. Part. Phys. 31 (2005) 873. 

[54] J. M. Allmond et al. Phys. Rev. C 78 (2008) 014302. 

[55] Coral M. Baglin, NDS 84 (1998) 717. 

[56] G. A. Lalazissis, S. Raman and P. Ring, At. Data Nucl. Data Tables, 71 (1999) 1. 

[57] R. Bijker, R. F. Casten, N. V. Zamfir, E. A. McCutchan, Phys. Rev. C 68 (2003) 064304. 

[58] M. A. Caprio, Phys. Rev. C69 (2004) 044307. 

[59] D. Bonatsos, D. Lenis, E. A. McCutchan, D. Petrelhs, I. Yigitoglu, Phys. Lett. B649 (2007) 
394. 

[60] S. G. Nilsson, Mat. Fys. Medd. K. Dan. Vid. Selsk.29, no. 16, (1955) 1. 

47 



[61] A. A. Raduta, D. S. Delion and N. Lo ludice, Nucl. Phys. A 564 (1993) 185. 
[62] G. Palma and U. Raff Am. J. Phys. 71 (2003) 956. 



48 



